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ABSTRACT 

A method for calibrating data obtained by the Urbana sodium lidar sys- 
tem is presented. Fii»t, an expression relating the number of photocounts 
originating from a specific altitude range to the sodium concentration is 
developed. This relation is then simplified by normalizing the sodium 
photocounts with photocounts originating from the Rayleigh region of the 
atmosphere. To evaluate the calibration expression, the laser linewidth 
must be known. Therefore, a method for measuring the laser linewidth using 
a labry-Perot interfere, -eter is given. The laser linewidth was found to be 
6 i 2.5 pm. Problems due to photomultiplier tube overloading are discussed. 
Finally, calibrated data is presented. The sodium column abundance exhibits 
something close to a sinusoidal variation throughout the year with the winter 
months showing an enhancement of a factor of 5 to 7 over the summer months. 
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1. INTRODUCTION 


1.1 iktru'vnl Li lav Cnnatpta 

The word Ildar is an acronym for Light l>et r- ‘Ion And Ranging. A* the 
name would imply, a lidar system is very simila. to a radar system operating 
at much greater frequencies. Like a radar system, a lidar system is composed 
of a transmitter and receiver. A lidar transmitter consists of a pulsed 
laser source sometimes equipped with benm collimating optics. The lidar 
receiver is comprised of a telescope, a filter to eliminate background light, 
a detector and signal recording equipment. 

I'or something to be detected by lidar methods, the object must interact 
with the light in such a manner that enough of the incident light is scat- 
tered back to the telescope that it can be detected above the system noise 
level, i'or atmospheric lidar, the scattering process is performed by parti- 
cles ami constituents which normally reside in the atmosphere. As a radar 
target's ef feet iveness is measured by a radar cross section, an atmospheric 
particle's ability to produce ar. echo is given by the particle's scattering 
cross section. The range from the lidar system to the scatterer is found by 
measuring the amount of time the light pulse takes to propagate to the scat- 
terer, be scattered and then return to the detection equipment. The range 
ot* the scatterer is given by '• ■ ~W.‘ where • is the speed of light and Lt 
is the time interval. The factor of 4 accounts for the two way propagation. 
It is assumed the scattering process occurred instantaneously, which is not 
always the case. An example in which the scattering process i r. not instan- 
taneous would he a resonant scatter from (T), at 10.4 im. The decay rate of 
this molecule for this transition is around ~ sec. 

This chapter will serve as an introduction by briefly reviewing lidar 
concepts and equipment as well ns a description of the sodium layer. Sec- 
tions l.J an<l l.S describe the historical progress of lidar systems and a 
summary of observed features of the layer respectively. Section 1.4 de- 
scribes the configuration of the present lidar system located near Urbana, 
Illinois. The chapter concludes with the objectives and scope of this study 
in Section l.f>. 

l.J .’•*/,?►' v>< r a(tt ''bser'at 

lidar, also known as laser radar or optical radar, was first used to 
study the upper atmosphere in the early l‘>ot)'s, soon after the birth of the 
laser. The first configurations were used to study t he atmospheric density 
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profile ns well as aerosol layers over a wide altitude range. The laser 
used with these systems was usually a Q-switched ruby laser. The operating 
frequency of the 1 nsmitter was not a critical factor in these observations, 
since the scattering mechanisms operate over a wide range of wavelength, A 
good review of lidar fundamentals and results concerning aerosol content and 
atmospheric density and composition can lie found in Kent and Wright [1970], 

With the advent of the tunable dye laser, scattering mechanisms which 
operated on an atomic or molecular resonance method could now be observed. 
Sodium, an atmospheric constituent which had been observed via airglow 
measurement techniques [Hunten, 1967], was the first species to be observed 
by resonance scattering [Fnu'nan et al . , 19f>9]. Several groups have become 
involved in the study of the sodium layer and hnve reported results concern- 
ing the nature and variation of the layer with the backing of an extensive 
data base \Pinonieh et a!., 1979; Gibson and Hand ford, 1971; Megie and 
R’.rant, 197*’]. 

As equipment improved and the integration time over which a profile was 
collected was reduced, reports concerning fine structure in the layer were 
reported [Kirohhoff tP'.d 'leen'sha, 1973; Hewlett et al . , 1978], A technique 
of processing the data using digital filters was developed which aided an 
observer in seeing the density perturbations in the vertical structure 
[.V/v.V? t et 1978] . There is considerable evidence which points to 

internal atmospheric gravity waves as being the origin of this structure. 

Other improvements in the lidar technique include a steerable system 
which can observe horizontal variations in the layer [1’homas et al., 1977], 
the study of the potassium layer \Megie et al . , 1978], daytime observation 
[Gilson .jnlPandfor<i, 1972] and observations of the vertical temperature 
profile [Planoni et al . , 1972; Gibson et al., 1979]. 

1.3 Po Hurt Loiter Pesoription and Reasons for P tu. In 

Atomic sodium is found over an altitude range of 75 to 105 km above the 
earth's surface. This altitude range corresponds to the upper mesosphere 
and lower thermosphere. The boundary between these two regions occurs at the 
mesopause. Inc peak of the layer is found somewhere between 83 km to 95 km 
with a width of from 5 to 12 km, The altitude of the peak can vary greatly 
over the spree of a few hours and is often bifurcated. Figure 1.1 shows 
data taken near Urban a on November 11-12, 1977 and shows a good deal of 
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variability. This data was filtered to remove spatial frequencies in the 
data higher than (3.5 km)”*. The first profile at 2120 CST shows a peak 
occurring at 83 km. The layer falls off very sharply on the bottomside and 
slowly decays on the topside. This small scale height, or distance over 
which the concentration falls by a factor of l/<?, is commonly observed and 
is much smaller than the atmospheric scale height. The bottom of the layer 
undulates up and down over the course of an evening having a period of about 
2 hours. At 235b, a second peak forms at 91 km, and in 30 minutes has become 
the dominant peak of the layer. These rapid variations and steep concentra- 
tion gradients point towards strong dynamic and chemical forces at play in 
this layer. 

Seasonal variations have been reported by several groups [Simoni ah at 
ah, 1979; Megie and Blamont, 1977; Gibson and Sandford , 1971]. They have 
reported increases in the column abundance during the winter by a factor of 
2 to 5 over those observed in the summer. Simonich et a!. [1979] postulated 

the variation in the ratio of increase was due to latitude, the stronger 
winter enhancements occurring at higher latitudes. 

The origin of the sodium in the atmosphere is thought to be from meteoric 
ablation. This is evidenced bv *-he number of other metallic elements which 
reside in this part of the atmosphere, each of which is commonly found in 
meteors. Other evidence of a meteoric origin is provided by observed en- 
hancements in the sodium column abundance during meteor showers, some as 
large at 200% greater than the monthly mean [Megie and Blamont, 1977]. 

There are three major reasons for wanting to study the atmospheric 
sodium layer, i ’ rst , the layer itself is important in the neutral and ion 
chemistry of the region. Although sodium is only n minor constituent at 
these altitudes, positive sodium ions are a significant species of positive 
ion and hence affect the charged particle concentration including electron 
densities. Sporadic-^ occurrences have been shown to be caused by enhance- 
ments of metallic ions, sodium included [Herrmann et at . , 1978]. Several 
authors have offered theories and models of the chemistry of the region with 
a recent one including cluster ion reactions [Richter and Reahrist , 1979], 

Secondly, the layer has proven to be an excellent tracer of dynamic 
effects in this region. Figure 1.1 shows many features which are thought to 
be due to dynamic effects, especially gravity wave propagation. The sodium 
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layer should provide information concerning these waves and I v they propa- 
gate, in the future. 

Finally the vertical structure of the atmospheric temperature can he 
measured in the sodium layer altitude range. There are two methods for 
making this measurement with lidar techniques. One wav is to step a narrow 
linewidth laser through 'he sodium line, collecting a photocount profile for 
each step. This results in an experimental lineshape at several altitudes 
for the sodium. Since t he lineshape of the sodium atoms is temperature 
dependent, the desired temperature information is deduced by finding the 
temperature which provides the host fit of a theoretical lineshape to the 
measured one. The second method involves placing an absorption cell in the 
received photons path. 'Hie absorption through the cell will be dependent 
on the atmospheric temperature. This measurement is important since tempera- 
ture information from this region is sparse. It is accessible for in nitu 
observation only by rockets. 

1 . -1 >’! p: ? " • rr' ’ •>. j'-' r.ir' U : > 

figure 1.2 shows a block diagram of the Urban a lidar system. The trans- 
mitter section is dominated by the dye laser. The dye laser used in the 
llrbana lidar system has been described by p ! [1978], 

However, a short synopsis of the laser will be presented here. This synopsis 
will also include new equipment added to the laser since the previous report. 

The tlve Rhodamine of! ,s used as the lasing material at a concentration 
of 5 * 10 S M. A solvent of equal amounts of cth;inol and distilled water is 
used mostly for safety purposes. A quenching agent, 1, 3, 5, 7 - cvclooctate- 
traene is used in the dye mixture at a concentration of 5 * 10 \M to aid the 
lasing transition's effectiveness. The dye is continuously circulated 
through a rod on the laser mainframe at a rate of 1.5 liters/min. 

The dye rod lies along the common focus of a bielliptical cavity. Two 
linear flash lamps having an arc length of 11 inches reside in the two inde- 
pendent foci. The flashlamps are the ablating type operating at a pressure 
of o - 9 Torr of nitrogen, fresh nitrogen is continuously passed through 
the lamps by bleeding nitrogen into one end of the flashlamps and pumping it 
out the other end with a vacuum pump. Ixcoss ablation products are removed 
from the system in this manner as well, being collected in a flask trap. 

The flashlamps are continuously water cooled to increase their operating 
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Figure 1.2 Block diagram of the Urbana lidar system. Double lines indicate light 
paths while single lines represent electrical connections. Angled 
lines in light path symbolize beamsplitters and mirrors. 
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life ami *o prohibit heat gradients from altering the index of refraction of 
the live rod, Since the flnshlamps create vibrations when they're triggered, 
care is taken to isolate mechanically the flash lamp setup from the optical 
rail where the dye rod and cavity optics rest, F.nergy for the flash is 
stored in a 5 id', 25 kV capacitor and is released by pro-ionizing the lamps 
with a high voltage pulse. The pre-ionizing pulse increases the conduction 
of the flnshlamps allowing the applied voltage from the capacitor to estab- 
lish a large electric field into the lamps then flash and release the 

stored energy. 

The dyi laser cavity must perform the selection of the operating 
frequency and sufficiently narrow the 1 inewidth of the output. This is 
necessary since the dye will 1 ase over a wavelength *ange of 5 nm. There- 
fore, the cavity has two frequency selective devices. At the rear of the 
laser, farthest away from the output end, a grating is placed oriented in 
the I. it trow position. This position implies the desired wavelength has 
its first order diffracted back along the cavity axis. Clearly, any other 
wavelength will be directed out of the cavity, thus providing the frequency 
selection. The current grating has 1800 lines/mm and is cut at a blare 
angle such that there is minimum loss into other orders. In front of the 
grating is located a flat 4 Of, reflective mirror oriented perpendicular to 
the optical axis. This mirror reduces the optical energies incident on the 
easily burned grating. These two optical elements together are called a 
mirror grating combination (MGC), The MOC provides a narrower 1 inewidth and 
higher reflectivity than the grating alone, aside from the reduced intensity 
on the grating. 

The next optical device on the rail is a Fabry- Perot etalon. An etalon 
is essentially two plane mirrors adjusted into high tolerance parallelism. 
This device further narrows the 1 inewidth of the laser and is the other 
frequency selective device. The two mirrors form an optical resonant cavity, 
having resonances at wavelengths where the round trip phase shift of a light 
beam within the cavity is a multiple of lv . . The separation of these resonant 
wavelengths is denoted as the free spectral range. This concept and other 
concepts associated with Fabry-Perot devices are discussed in Chapter 5. 

The free spectral range of the etalon presently in use is 0 . (> nm. By 
changing the angle of incidence at which the beam passes through the etalon, 
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the effective length of the resonant cavity is changed and the laser fre- 
quency is altered. Fine tuning of the laser is effected in this way. On 
the other side of the dye rod, the output mirror is put into position. The 
current configuration uses a spherical mirror with a radius of curvature of 
A m and a reflectivity of 55%. A 45-degree folding mirror directs the beam 
upwards into the atmosphere and completes the rail optics. 

Tuning is monitored by a grating spectrometer and a sodium resonance 
cell. The coarse tuning is performed by using the grating spe .trometer with 
a sodium arc lamp providing a D j and line for reference. Unfortunately, 
the spectrometer does not have enough resolution to tune the laser output 
onto the 3 pm wide sodium line. This fine tuning is done by passing a por- 
tion of the laser beam through a heated cell filled with sodium vapor. When 
the laser is at the correct frequency, scatter can be observed as the beam 
passes through the cell. By maximizing scatter, good tuning can be assured. 
A picture of the laser mainframe is shown in Figure 1.3. 

The receiving system begins with the 1.22-m Fresnel lens telescope. 

This instrument is described by Rowlett and GdViiner [1979] in detail. A 
brief description of it follows for completeness. The telescope consists 
of the 1.22-m Fresnel lens, a folding mirror to change the optical axis from 
vertical to horizontal, a field stop iris, a collimating lens, a filter, a 
focusing lens and finally, a photomultiplier tube (POT). The collimating 
lens is necessary to allow proper performance of the interference filter. 

The POT is a high gain, low noise unit. The system noise level is dominated 
instead by background light collected by the telescope. 

The output from the POT goes to the laser radar interface [Kin ter, 
1977]. In the interface is a high speed counter, continually being incre- 
mented by pulses from the POT. The interface begins timing once it has been 
notified by the laser pulse detector that a shot has begun. After a delay 
period when the returns begin from an altitude of interest, the counter 
is periodically sampled. F.ach sample is stored in a quick access memory 
location. Once the shot is completed, the data is transferred to a Z-80 
microcomputer [ I'eitelbaum and Sechriet 3 1979] where the difference between 
samples is computed and stored. The sampling interval corresponds to a 
height interval from which the counted photons originated. Therefore, each 
value stored in the computer represents photons from that range interval or 
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range bin. In this way, an altitude vs. photocounts profile can be obtained. 
Very few photons originating from the sodium layer are obtained by one shot. 

To construct n good, low shot noise profile, range bin values over several 
shots must he added together or integrated. Usually anywhere between 50 and 
300 shots are required for each profile spanning an interval of 5 to 12 
minutes of integration time. 

Once enough photons for a profile are collected, the data values are 
stored on cassette tape. After the evenings observat ion is over, the data are 
re- rent! into the computer and transferred over the phone line via a modem to 
the University of Illinois Cyber 175. At a later date, the data are processed 
in final form. A table of the parameters of the Urbana lidar system is pre- 
sented in Table 1.1. 

1,5 I’/vV.'fi'iv’p Th;p i'tu.fy 

The objective of this study is to make the necessary measurements and 
calculations necessary to convert the sodium photocount data into sodium 
concentration vs. altitude profiles. The concentration value are impor- 
tant factors in the validation of dynamic and chemical models of this region. 
Calibration also allows profiles to be compared within an evening, from day 
to day or on a seasonal basis. 

The major obstacle preventing the calibration of this system in the 
past was an unknown value of the laser linowidth. Therefore, a major part 
of this work was involved with this measurement. Other factors, such as PNfT 
overloading, caused previous methods to he outdated. Finally, the data had 
to he processed to obtain calibrated data. 

Chapter 2 presents the principles and concepts behind the calibration 
technique. Chapter 3 provides information concerning the linewldth measure- 
ment and the theory behind it. chapter 4 summarizes the TMI’ overloading 
problem while experimental results and error analysis are given in Chapter 
5. Conclusions and recommendations for further study comprise Chapter b. 


Table 1.1 


Values of t ho current Urban a liilar system 
parameters . 


Parameter 


Value 


Wave leapt h 

Transmi 1 1 ed Ramliv i ilt h 
laser Beam Piverpenee (half an pi el 
l aser Pulse I nerpy 
Pulse 1'urat ion 
Receiver onto Time 
Receiver Area 
Receiver Bandwidth 
Receiver I ffieieney 
Re ce i ve r Be am I ' i ve rpen c e (half an pi e ) 
Separation of Receiver and Transmitter 
Repetition Rate 


S8R nm 
f> pm 
*■ 1 mrad 
SO- 100 mJ 
Uts FWIIM 
Sus 

1.1" m“ 

S nm 

o.os 

S mrad 
10 m 
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2. Lll)AR PRINCIPLES AND CALIBRATION TECHNIQUE 

2.1 Introduction 

The previous chapter gave a brief overview of lidar system operation. 
This chapter analyzes the process by which signal photons are received and 
counted, in enough detail to he able to quantify the process. A relationship 
between the observed photocounts and the concentration of the scatterer pro- 
ducing the signal photons is developed. The material is presented more for 
completeness than for originality, and the development closely parallels the 
one presented hv Richter and Sechriat [1978]. 

Section 2.2 discusses how the major loss mechanisms operate on a light 
beam propagating through the atmosphere and its relationship to the received 
signal. Section 2.3 develops the lidar equation, relating photocounts to 
sodium concentration. The effective scattering cross section and the trans- 
mission through the sodium layer are discussed in Sections 2,4 and 2.5, 
respectively. Hath of those terms are used in the calibration calculation. 

A method of calibration which simplifies the expression derived in Section 
2.3 by normalizing the sodium returns to returns from the Rayleigh scattering 
region of the atmosphere is the topic of Section 2.6. The chapter concludes 
with a description of the Rayleigh scattering cross section, a quantity used 
in the normalization process, in Section 2.7. 

2.2 Liaht Bcxrm Extinction in the Atmosphere 

There are two mechanisms which operate on a light beam to cause it to 
be attenuated as it propagates through the atmosphere. First, the light can 
be absorbed by particles with absorbed energy being transferred into some 
non-optical foim of energy. Secondly, the light can be scattered from the 
beam into another direction. The expression for the intensity of a beam 
after it has propagated through a media of length *, is given by Bouger’s Law 
which is: 


I - J 0 expl-J^ ^JeUs] 


( 2 . 1 ) 


where 8 is the extinction coefficient over the path of propagation. The 

C X 

extinction coefficient is made up of two parts, the volume total scattering 
coefficient 8 p , and the absorption coefficient 8^, in the following manner: 


ex 


' 8C 




( 2 . 2 ) 
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The term ^ is particularly important for lidar systems since it is the 

source of the signal pr tons collected by the system. The units for fl , 

h and p , are in inverse length, 
n » it 

There are four basic types of scattering processes which can contribute 
to e . They are Rayleigh, Mie, resonance and Raman scattering. Raman 
scattering is characterized by a frequency shift in some of the returns, 
ibis effect stems from a change in energy levels in an excited molecule and 
requires the quantum theory to describe. This type of scattering is not 
important for the Urban a lidar, and so will not be dealt with hereafter. 
Rayleigh and Mie scattering are caused by the presence of particles in the 
media o ** propagat ion which form inhomogenci t ies for the propagating fields. 
Rayleigh scattering theory is a special case of Mie theory, applicable when 
the radius of the scattering particle satisfies the following relationship: 

>• <_ . 03V ( J . 3 ) 

Ibis condition is satisfied In* gaseous atmospheric molecules which comprise 
the air. The Mie theory is applicable to particles such as volcanic ami 
meteoric dust as well ns water droplets in clouds. Rayleigh scattering is 
characterized by equal amounts of scattering into the forward and rear hemi- 
sphere of a particle. I he forward direction is determined by the incident 
he.uiis direction. The Rayleigh scattering at right angles to the beam of 
incidence is strongly polari.ed. Mie scattering exhibits strong scattering 
into the direction of propagation. The larger the ratio of radius to wave- 
length, the stronger the forward scattering will be. 

Finally, resonance scattering involves the absorption of n photon by a 
molecule nr atom into an electronic transit ira and its subsequent spontaneous 
emission, I'he emitted photon is nt or near the original frequency and 
travels with equal probability in all directions. This is the mechanism h\ 
which the sodium layer is observed. 

I'o he able to quantify further the effectiveness of a scatterer, some 
assumptions must first he made, for atmospheric scattering, the particles 
are separated adequately and undergo random motion so that the contributions 
to the total scattered energy by many particles have no phase differences. 
This allows the total intensity to be calculated by adding the intensity 
from each scattering particle. This process is known as independent scat- 
tering. Vnother assumption made due to the low density of air is that of 
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single scattering. Single scattering implies that a photon is scattered only 
once. If a photon is scattered two or more times then multiple scattering 
occurs. Naturally, single scattering results in a much more straightforward 
manner of describing the scattering process. 

Several quantities exist which aid in quantitatively describing the 
effectiveness of a scatterer. The angular scattering cross section, o(9) is 
defined as the area of impinging flux which is scattered into a unit solid 
angle in an observation angle, fl. This geometry is diagrammed in Figure 2.1. 
By this definition, the number of photons scattered from a photon flux ♦ 
into direction fl is given by: 

.V * i> ofR) photons per steradian (2.4) 

The total scattering cross section is another quantity used in scattering 
theory and is defined as the total area of flux scattered in all directions 
by the particle. This can be found by integrating rr(e) over 4* steradians. 

Two more commonly used terms are the volume angular scattering coeffi- 
cient and volume total scattering coefficient. The volume angular scattering 
coefficient B ( o ) is defined as the area of incident flux which is scattered 
in a direction 0 per unit solid angle by a unit volume of scatterers. Under 
the assumption of independent scatter, this quantity would just be the 
angular cross section of one scatterer multiplied by the total number of 
scatterers in the unit volume. This can be expressed as: 

B(o) * o(fl Ip (2.5) 

where p is the number density of scatterers. The volume total scattering 
coefficient is the term used in the beginning of this section in Bouger's 
Law. Tt is defined as the area of incident flux which is scattered in all 
directions by a unit volume of scatterers. A value for p can be found by 
either multiplying o by the numb-'-" density of scatterers or by integrating 
B(9) over 4n steradians. 

Hopefully, these definitions will provide the needed groundwork to 
understand better the development of the lidar equation in the next .section. 
Sections 2.4 and 2.? will discuss in more detail values for .a resonance and 
Rayleigh cross section, respectively. Several good references exist to de- 
scribe scattering in the atmosphere. McCartney [1976] is a thorough and 
highly readable text on the subject while Van De tiulst [1957] is a rigorous 
and often cited resource. 
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2 , 1 !i lor b'qifTt ?' rm 

The lidar equation relates the received signal obtained from a scat- 
tering region to the concentration of the scatterer as a function of several 
system parameters. There are numerous ways to develop this expression. A 
very rigorous method is presented hy Ma.wurva (1977J. Another derivation, 
given by r' *h *rr rv:.‘'f 'hri8* [ 1 0”* k ] presents a more intuitive development 
and gives a result which is applicable to the Urhana lidar. Thnt development 
will be closely followed in this section. 

First, consider the geometry of the system as depicted in Figure 2.2. 
This system operates in the monostatic mode in which the distance of separa- 
tion between the transmitter and receiver Is negligible compared to the 
distance at which the photocounts originate. Due to the method the counter 
on the photomultiplier tube is sampled, the height is divided into intervals 
of altitude called range bins. The range bln at the bottom of the sodium 
Inver will he denoted as range bin 1. Tito altitude at the bottom of range 
bin 1 is , with the center of the range bin occunin? at an altitude of 

n * A r.'2. Here, represents the thickness of the ringe bin and is related 

* j 

to the counter sampling time bv A;: * ^ .’Af where .» is the speed of light, 
rhe notation is easily extended to range bin t for the general case. 

Cons i do i now a photon flux impinging on the bottom of range bin i. The 
flux expression is a function of frequency v and represents the number of 
photons per unit area within a frequency interval v to v ♦ <iv. The flux can 
he expressed as: 


t( v 


fl.) 


V ;7(V) 

“ A. in) 7 


( 2 . 6 ) 


where Y, is the total number of transmitted photons, < 7 (v) is the distribution 
of the transmitted light over frequency, T is the transmission of the beam 
through the atmosphere to height it., and A^z^) is the area illuminated by 
the laser beam at this altitude. The distribution function p(v) is normalized 
and defined bv: 


r 

j n 


r(vl :v = 1 


(2.7) 


The term p(vt is the laser lineshape in this case and when evaluated at v, 
represents the fraction of the total nuiiiher of photons found in the frequency 
interval of v to v + :v. 



lUajtran shouin)’ t!u t ransnu 1 1 ed lijjht's eeomet rv , j no 1 ikI i n> 
several ranr.e bins at different heights. 
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To write equation (2.1), it was assumed the laser pulse can approxi- 
mated by an impulse function. That is, all of the photons in the laser pulse 
are assumed to be at the same altitude and impinge on the bottomside of range 
bin i together. In reality, this is not the case. A more rigorous develop- 
ment will show the total intensity received will be proportional to the two- 
fold convolution of the scattering medium's concentration with height, the 
pulse shape of the laser in time, and the receiver gating function. The 
receiver gating function describes how the photons are sampled. For the 
(Jrbana lid' , this is a square pulse having a value of 1 during the time 
between samples and zero otherwise. As is the case with a convolution, the 
gating function and laser pulse shape will tend to obscure resolution in the 
concentration profile by smearing out the profile's fine structure. This 
process of obscuring fine structure as well as a development of the two-fold 
convolution can be found in Richter and Seohriat [1978]. 

Fortunately, through the use of digital signal processing means are 
available to reduce the effects of this convolution. Rowlett and Gardner 
[1979] discuss a technique for filtering the received signal with digital 
deconvolution filters. This work has proven quite successful in recovering 
the sodium concenti at ion profile. Since the effects of a finite width pulse 
and gating function are mostly removed by this process, they will be neg- 
lected in this derivation. 

For the resonance scattering process an absorption and spontaneous 
emission are required by each scattering atom. If the atom is collisionally 
deactivated, it will not add to the scattering process. Also, if the decay 
rate of the atom is too slow the range calculation will become more involved. 
For the case of the sodium layer, neither of these processes are important. 
The reason for this is that the decay time for the D ^ transition is 16 nsec 
which is much smaller than the sampling time of 3 usee. This time is also 
much smaller than the average period between collisions at these altitudes 
which is between 1 to 100 usee. Therefore, there is a negligible probability 
of the atoms either taking so long to decay that the resolution of the system 
is impaired or that the atoms will be collisionally deactivated. 

An expression can now be developed tc predict the number of photons col- 
lected by the receiving system. Figure 2.3 outlines the associated geometry 
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for thin collection process. Recalling the definition of the argular scat- 
tering cross section, the number of photons per frequency interval scattered 
by one sodium atom into a unit solid angle is given hv: 

v , t (0 * n(i> ,v (2.8) 


The e dependence of .V can be removed by recalling that resonance scattering 
is isotropic. Therefore, rr(0,v,;:) can be replaced by the total scattering 
cross section divided by U steradians to obtain a unit solid angle basis, 
liquation (.’.H) can now be rewritten as: 


l v 


Q ( v , r . ) 
•In 


(2.9) 


The number of sodium atoms in the differential volume shown in Figure 
2. 3 wi 1 1 he: 


Jo) 0(a) :T(:: I 


1 ( •' ) ,! ( in)::: 


( 2 . 10 ) 


where , is the number density of sodium atoms at height r.. It is assumed 
here that the laser divergence angle is small enough that the volume of the 
space which comprises the range bin car lie approximated as a cylinder. The 
number of received photons can now be calculated as follows: 


Tv ,; 


(o) 


V flt (v,al a (a) "n 


( 2 . 11 ) 


where is the solid angle subtended by the receiver aperture from a scat- 
tering atom, is the atmospheric trnnsmi ss ion bach through the atmosphere to 
the receiver, and n is the receiver efficiency. The solid angle subtended 
by the receiver is given by: 


( 2 . 12 ) 


The proper substitutions may now be made into equation (2.11) to obtain: 

. \\ .7(v) A 

Y„( (.'(••) x ;) * (_£) x 7' x n (2.13) 

r.ineelling out .,(;:) and rearranging yields: 



"i * 




(v 


(2.14) 


Integrating each side over frequency in (2.14) removes the frequency depend- 
ence of and results in: 


21 


j) * - 


N.A T n , ^ 

t r £l$l 
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• 0 


$(v) a{\>,z)d\> dz 


(2.15) 


The term in brackets is defined as the effective cross section, a 


V.' (a) 


a(v) a{v,z)dv 

Jo ‘ 


(2.1b) 


The z dependence of cs rr is due to temperature variations over altitude. For 
large enough linewidths, the temperature dependence is negligible and hence, 
so is the dependence on altitude. This approximation as well as other fac- 
tors pertaining to the calculation of ^ are discussed in Section 2.4. 

Finally, equation (2.15) can be integrated over the altitude interval 
of range bin i to obtain the number of photons one would expect ♦‘o count 
from that range bin. Performing this integration produces the following 
expression : 


.V. 

1 


.\\a : 
( v 


4n 




'-.1. 

s" 


(2.17) 


Since the interval of integration is small, p(n) can he approximated by the 
average concentration in range bin i , p. and s can be approximated by the 
average altitude of range bin i , namely (z^. + Aa/2). Sodium layers have 
shown scale heights as small ns 2 km [Gibson and, Gandfoni, 1971] which is 
still fairly large compared to the range bin thickness of 450 m. The range 
bin thickness is also only a small percentage of its altitude. Therefore, 
these approximations will result in very small errors. Upon making those 
approximations, the integration becomes trivial allowing equation (2.17) to 
be expressed as: 


.V. = 

t 


•VV-'-n p . 


Ljl.* 


4ir(a^. + A z/2Y 


(2.18) 

A value for 


Recall that T is the transmission to the bottom of the layer, 
it will be developed in Section 2.5. 

2.4 J v 'Tvf nv Font tcrnm dross deal ion 

The effective scattering cross section as defined in equation (2.1b) 
relates the scattering response to the frequency dependence of the atom as a 
scatterer combined with the frequency distribution of the light being used 
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for excitation, fissent i al lv , equation (2.16) convolves the scattering cross 
section of the atom with the laser line shape. 

The sodium resonance line shape is shown in Figure 2.4. Here, both the 
and P, lines are presented. Closer examination of the P 2 line shape shows 
that it is composed of two groups of three hyperfine lines. These groups are 


denoted as the and line for sodium. 


The P_ line is located -2.14 
2 a j 


m from the line center position while the P line is 3.56 m away. Each 
of these lines is Doppler broadened and will have a linewidth given by: 


Av, 


2 iZF 1 n 2 


(2.19) 


where : 


' p c o y 

* universal gas constant * 8.314 J/mole 
« speed of 1 i ght 

v (1 * center frequency of the broadened line 

7 » temperature 

y * molecular weight * ('.23 hg/mole. 

Doppler broadening is described by a Gaussian line shape. Therefore, 
the scattering cross section for the entire P-, line can be expressed as the 
sum of two properly weighted Gauss ians. This weighting cal ail at ion was 
performed by h’i.’htrr' 7>: ! [19"8] and results in a scattering cross 

sect ion of: 


cr(v) * a 
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(2.20) 


Here o (1 is the peak scattering cross section for the combined line and has 
a value of 


1 . bt> x HI 
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Av. 


m 


( 2 . 21 ) 


The terms and v,^ represent the center frequencies of the P 0 ^ and P^ 
lines, respectively. The variance of the two Gaussian expressions o^, is 
defined as: 


a,, 


Av,," 

8 ln2 


( 2 . 22 ) 


The lineshape of the laser can also be represented by a Gaussian of the 


form: 





24 


.»lv) * — PXp - 

a, - 
Av ; 

where ^ -y— y ■ variance of the laser line 

v fJ » laser center frequency. 

*v. * laser 1 inewidth. 



(2.23) 


Substituting equations (2.20) and (2.23) into equation (2.16) yields: 


0 

t 


> /*/* 




exp 


-> 



20 

33 exp 





12 

♦ 33 exp 



The 

_oo c; 

That 
by :i 


lowor limit of intonation on equation (2.24) can he extended from 0 to 
inco there is a negligible contribution to the integral from this interval, 
transforms the above integral into a convolution which has been solved 
pi i ;'c*'hv*rt [1978]. This new expression for a at y is then: 
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— 1 } (2.21 

♦ Av“ JJ 


Therefore, by knowing the temperature of the sodium as well as the linewidth 
and center frequency of the laser, the effective cross section may be calcu- 
lated. Figure 2.5 shows how depends on the linewidth for two values of 

the laser center frequency for sodium at a temperature of 200 K. The dip in 
11 . f° r lasers tuned to line center for small linewidths corresponds to the 
dip in the line at line center as observed in Figure 2.4. Clearly for 
small linewidths, tuning becomes very critical to obtain maximum values of 


The previous section made 
altitude. The validity of that 


the assumption that a ^ did not vary over 

^ i 

assumption may now be investigated. Figure 
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2.6 shows the temperature dependence of a a ^j> for linewidths of 0.1, 1 and 10 
pm for lasers tuned to line center. A glance back at Figure 2.4 shows that 
the cross section figure is much more sensitive to temperature here than 
anywhere else. At midlatitudes, temperatures can vary from 170 K to 220 K 
over the altitudes where sodium is found [7/.S. Standard Atmosphere , 1976]. 
For a large linewidth such as 10 pm, little if any variation can be observed 
over this temperature interval. liven for linewidths as small as 0.1 pm, a 
variation of around . + 10\ from the 200 K figure is evident. Therefore, it 
seems that assuming to be constant over altitude is quite reasonable. 

One lost approximation can be made to equation (2.25) for large line- 
widths. Av. dominates Av n allowing Av. to be neglected in the linewidth sum 
terms. This results in: 

-6 


1.66 x l'; 
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The llrbana dye laser has a linewidth of 6.5 pm which equates to 5.62 x 10 Hz 
while the Doppler broadened line has a width at 200 K of 1.07 x 10^ Hz. 
T:quation (2.26) yields a value less than 2% in error for o ^ compared with 
the more accurate expression in equation (2.25). This implies that equation 
(2.26) can be used with good accuracy for the present dye laser configuration. 
2.5 Transmission Losses Through the Sodium Layer 

In Section 2.3 the entire losses incurred by the beam as it propagated 
from the ground to range bin i were lumped into one term, T, Consider 
breaking this transmission coefficient into two terms. One term denoted T^ 
represents losses through the atmosphere to the bottom of the sodium layer. 

The second term will account for the transmission through the sodium layer 
only and will be denoted as The total transmission term can now be 

represented by: 


A 

«*t 


Wrt 


(2.27) 


Before computing a value for , first consider a beam's transmission 
through an arbitrary range bin, where both the sodium concentration p and 
the cross section o(v) are constant throughout the range bin thickness. 
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Tor an impinging beam with a lineshape of <?(v) , the number of photons 
absorbed by the layer in a frequency interval of v to \> ♦ can be expressed 

cs : 


■Vj(v) * A’ 0 j7 (v) [ 1 - T’(v)] ■ N q p(v)[l - exp(-c(v) pAa)] (2.28) 


where is the number of absorbed photons per frequency interval and A'^ is 
the total number of photons at the bottom of the range bin of thickness Aa . 
The absorption is considered here due to its expression being easier to work 
with in succeeding steps. To obtain the total number of absorbed photons 
over all frequency, each side of equation (2.28) must be integrated over all 
frequency. This is done by first assuming that both p(v) and a(v) can be 
expressed as (Russians. The preceding section modeled a(v) as a sum of two 
(Russians, however a single, equivalent Gaussian should provide accurate 
results for lineshapes with a linewidth equal to or greater than the com- 

ft* 

plex width of 7> pm. Expressions for p(\0 and afv) will then be given by: 


a ( \) ) * 





*(v) 
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-U) 


(2.29) 


where: » linewidth normalisation constant 

a n m maximum cross section value * + a., 

0 On Or 

u> * 2>ln2 (v v () )/Av,, 
a * Av,/Av„. 

Ij . ' 


This notation is somewhat different than that previously used, hut will aid 
in subsequent steps. 

Under these assumptions, an absorption coefficient, A can now be calcu- 
lated in the following manner: 


.V a ( M )./w 

A . _£ ( 2 . 30 ) 

•Y 0 p(u>) iw 


The usual assumption concerning limit extension has been invoked here. This 
method of obtaining A is chosen to allow the integration to be carried out 
with respect to w without having to go through the formality of a change in 
variables. Substituting (2.28) and (2.29) into equation (2.30) yields: 
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j~ { 1 - exp[-<r 0 #’ u) pA*] } «’^ a ^ <fa 
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4' Ma] du 


(2.31) 


This quotient has a series expansion [Mitchell and Tmaneky, 19~1 ] 
given by: 
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0 0 pA* 
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n ,/T7 


(2.32) 




Assuming values of a n * 1.66 x 1 ()”*Vav„ m 2 , p * 5 x 10^ m Aa ■ 450 m, 

q ” * q 

Av. - 5.62 X 1 O' |I 2 (AX , - 6.5 pm), Av M « 1.07 x 10* Hz (T - 200 K) yields 

-4 ” * -0 

values of 6.53 x 10 far the first term and 8,11 x 10 for the second term 
in the series. '’learly the first term will dominate the series allowing the 
higher order terms to be dropped from equation (2.32) without much loss of 
accuracy. A value fo: the transmission through one range bin can now be 
expressed as; 


1 - A * l 


o Q pAn 


[1 * (Av^/AvJ 2 V J 


Combining the effects of all bins below bin i , the transmission of the beam 
from the bottom of the layer to bin i will be given by: 
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(1 «• (Av r /AVp) 2 ] 5 


(2.34) 


hquation (2.34) represents the transmission of the laser beam upwards 
through the sodium layer. The light that travels downwards does not origi- 
nate from the laser, but from the atoms themselves. Therefore, the distri- 
bution in frequency of the downwards travelling light will be different, 
causing a different transmission factor to be in effect. Assuming the 
distribution of the downwards travelling light can be approximated by a 
f.aussian of linewidth Av n , equation (2.34) will still be in effect after 

replacing Av, by Av.,. The downwards transmission through the sodium layer 
* / * 

can be expressed as: 

i- 1 a . p . A r, 

7 . * :i (1 - ) (2.35) 
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The upwards transmission ns given in equation (2.34) will now he denoted as 
T . Applying this terminology to equation (2.18) viclds: 
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(2.36) 


One might wonder if the detailed expressions for and T. nr.* worth 
the effort to calculate. To determine their effects on the concentration 
calculation, a simulation was carried out. The count values were assumed to 
have a C.aussian distribution centered at 90 km given by: 


’ 0 *' 


(a 
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(2.3-) 


Tiie terms .V and , 7 " as well as a lumped value for all the constant terms in 

1 9-3 

eq .ait ion (2.36) w'*ro chosen to yield a concentration of 5 * 10" m ' at 90 km 

1 3 2 

and a total column abundance of 5 « 10 m . These assumptions resulted in 
the concentration having the following representation: 


r . = 6.2 * 10 


, (r.. * .225)' 


v a. - 90)' 

Ton 


(2.38) 


where a. is given in kilometers and 0 • in m ' . 

The results of the simulation are shown in figure 2 .". The figure shows 
the full correction, for both upward and downward propagation, accounts for 
;ui inert se in the peak concent rat ion of 10*. Most of that contribution 
comes from the downward propagation correction as revealed by the figure. 
Therefore, it seems worthwhile to include the sodium layer transmission terms 
in the expression for the sodium concentration. 

2.6 .Vcmt/ fr-tr * ■*»; Mi'tk' : of ibmt ion 

In the expression for the counts from range bin z as given in equation 
(2,36), there are some factors which have yet to be discussed. Two of these 
factors in particular may prove difficult to measure or calculate. Those 
factors arc the atmospheric transmission, , and the total number of pho- 
tons, V,. At present, a nightly value of ”, does not exist for the llrbana 
lidar. This term should ''ary considerably for the different nights of obser- 
vation due to different levels of haze, aerosols, clouds and other components 
which cause attenuation of the laser beam. It straightforward to measure 
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the total number of photons. However, previous measurements have shown this 
figure to change on a shot to shot basis. Fortunately, a method exists which 
can eliminate both of these terns. 

The elimination of difficult to evaluate terms is done by normalizing 
the returns from the sodium layer with returns from the Rayleigh region of 
the atmosphere. The returns from the Rayleigh region have been well studied 
and car. be accurately predicted. To see how this normalization is carried 
out, an expression for the number of photocounts from the Rayleigh region 
must first be derived. 

Consider a photon flux impinging on a range bin located in the Rayleigh 
scattering region of the atmosphere. The thickness of this bin is As and 
the cer.:er of the bin is located at an altitude r,^. Rayleigh scattering is 
not nearly as sensitive to rrequency as resonance scattering. Therefore, the 
frequency dependence of the 1 user lineshape need not be taken into account. 
Under this assumption the photon flux can he expressed as: 

■V, 

♦ (a) r - 7 - -r? (2.39) 


where T h is the transmission to the bottom of the Rayleigh bin. The number 
of photons scattered by one air molecule per unit solid angle towards the 
telescope will be: 
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(2.40) 


The observation angle a is tt radians for this case since 6 is measured from 
the initial direction of propagation. Backscattering is the term applied to 
this geometry. The number of scatterers in a differential volume will be: 

>■;, = P r A, (2.41) 


These equations are combined in a manner similar to equation (2.11) to obtain 
a value for the number of photons counted by the system from the differential 
volume : 

&' r (a) = n p (z) A 3,00 ft(s) '” R n (2.42) 

where n and ; have the same definitions as given in Section 2.3. Making the 
proper substitutions , integrating over the altitude of the range bin, and 
making similar approximations to those made in Section 2.3 to evaluate the 
integral results in: 
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It was assumed that o,,( n ) was constant over the range of integration. Next, 
dividing equation (2.3(0 by equation (2.43) and solving for p ? . yields: 
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(2.44) 


It was shown by Rii'htrr m l >\' 'hrist (1978] that for values of stj, of 30 km 
or greater scattering and absorption losses to the beam from this altitude 
to the bottom of the sodium layer were about 11> in the absence of aerosols. 
Since significant amounts of aerosols rarely occur at this altitude, and 
", may be equated without much loss of accuracy. It is common in sodium 
lidar work to incorporate the factor of 4n into the Rayleigh angular scat- 
tering cross section, liquation (2.44) may now be rewritten as: 
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(2.45) 
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This is the final form of the calibration equation. The next section dis- 
cusses how a value for rr.i-n) is obtained. Chapter A will discuss how an 
altitude, r r , is chosen which will provide pood results. 

J,” luiUit'fah ;\\it ft'tr na < ’rvns 

lo derive an expression for the an pillar scattering cross section, first 
consider the volume total coefficient ft. Recall that ft represents the amount 
of energy scattered in all directions by a unit volume of scatterers. An 
expression for it is piven as [A , 197b]: 
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(2.40) 


where »: and e. are the index of refraction for the scattering media and p 
is a constant determined by the composition of the scatterer. The term in 
brackets is a depolarization factor due to molecular anisotropy with p 
having a value of U.035 for air [Pirmufcvf, 1HS?]. 

Next, the concept of phase function needs to he introduced. The phase 
function is defined as the ratio of energy scattered into a direction per 


unit solid angle to the average energy s- attered in all directions per unit 
solid angle. Under this definition the phase function r(6) must satisfy the 
following condition: 


!_ 

4ir 


4 TT 
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r(e) <iu - 1 


(2.4?) 


where L is a differential solid angle. Anisotropy also affects the phase 
function. An anisotropy corrected F( O) for air is given by McCartney (1976): 

F(n) « 0.7029 (1 ♦ 0.9324 cos 2 9) (2.48) 

This phase function is for unpolarized light, fiach of the components inside 
the brackets is for a particular component of polarized light. The 1 term 
is for light polarized in a direction perpendicular to the observation plane 
while the cos “a term is for light polarized parallel to the observation 
plane. The observation plane is the plane which contains the line indicating 
the direction of travel for the incident light and the line in the direction 
of observation. For baekscattering, each component contributes fairly 
equally, so that details concerning the polarization are not necessary for 
an accurate value of the phase function. 

With the aid of the phase function, the angular volume scattering 
coefficient may be derived. Since the phase function represents the ratio 
of the energy scattered into one direction per steradian to the average 
energy scattered in all directions per steradian, the angular volume scat- 
tering coefficient may be found by multiplying P(e) by the average energy 
scattered in all directions per steradian. This quantity is simply P divided 
by 4 it steradians. Performing this multiplication results in: 
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0.7629 (1 + 0.9324 cos~e) 


(2.49) 


Recall that in Section 2.2, the angular scattering cross section could be 
found by dividing p(,0) by p.,. Remember also that the term a,,(ir) as used in 
Section 2.6 included a factor of 4* . Performing these manipulations results 
in : 
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0 . ”629 (1 + 0.9324 cos 2 0) (2.50) 
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!\ , wui?rf [15157] has presented good tables by which (n • 1) might he calcu- 

i i 

lnted. The (u“ - 1)“ term mny be expressed in terms of (n - 1) since n is 
so close to 1 for air. This approximation is given as: 


(n 2 - l) 2 s 4 in - l) 2 (2.51) 

Since o p (fl) is only used for air and a * it , the depolarization factor and 
phase function may also he evaluated to combine with equation (2.511 to 
transform equation (2.50) into: 


32i» (*: - 1)“ , „ . 16n 3 (« - 1) 

a,.(°) S— t l- 53 ' S — j- L 

il , »-»*♦ s “ 


3f\ X 


P K * 


(2.52) 


Another way also exists to calculate a value for o^,(0) . MaCtatchey and 
p'Aaxti [1978] gave an expression for 8 with corrections for frequency pres- 
sure and temperature as: 
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1013' 4.0117 


(2.53) 


where 7 is in Kelvin, P is in millibars and X in centimeters. The units of 
8 arc inverse meters, The dependence on wavelength of instead of X* 

is line to (>: - 1) having a slight dependence on wavelength. Following the 
procedure used earlier in the section, o p (ir) can he calculated by 8 in the 
following manner: 
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1.4379 


(2.54) 


Substituting equation (2.53) into the expression results in: 
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1.41 x 10' 


273 V 
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,4.0117 


(2.55) 


Figure 2.8 illustrates the change in it._,(it) over an altitude range of 25 to 50 
km. Values for P. and p„ were taken from the IKS. Standard Atrvosphcw 

n 

[197(1 ]. As can he seen, a v ( tt ) varies only S*o over the plotted altitude range. 
Deviations in the standard atmosphere parameters from the actual atmosphere 
arc too great to make the variation of o^(tt) with altitude significant. 
Therefore, Op(tt) as evaluated at 35 km will be the value used in this study: 

a.,(ir) = 4.75 * 10" 31 nf 2 (2.56) 

Tliis value is between that reported by Bidhtcr and dedhvi-at [1978] and the 
value reported by [1976], and can be assumed to incorporate the latest 

evidence available in the literature. 
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3.1 Introduction 

The previous chapter developed an expression for the effective scat- 
tering cross section, which is reproduced here: 
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(3.1) 


The cross section is a function of only two parameters which can vary, the 
laser center frequency and the laser linewidth Av^. Both of these quan- 
tities are in terms of frequency units but are related to wavelength units 

■> 

by » i’/Vp and AX^ * The wavelength units convention will be 

used throughout this chapter. 

The Urban a liilar system uses a sodium resonance cell as described in 
Fichtcr ind do 'hriot [1078] to monitor the position of the center frequency. 
To briefly describe this, the laser beam is passed through a cell which 
contains amounts of sodium vapor. The cell is heated to a temperature of 
around ISO C to provide a large enough partial pressure to ensure good 
resonance scattering of the beam when it is at the proper wavelength. When 
the laser is "off line", the cross section of equation (3.1) is very low 
and the beam passes through the cell unattenuated. When approaches line 
center of the sodium complex, the cross section increases and one is able 

A* 

to observe scattered light in the cell. Once the laser has been tuned "on 
line", it is monitored frequently through the night to ensure that it stays 
"on line". Drifting due to thermal fluctuations and other factors does 
occur, so that adjustment is required from time to time. This study will 
assume this method is sufficient to assure the laser frequency is at line 
center. Errors to this approximation will be discussed in Chapter 5. 

The second quantity that could vary in equation (3.1) was the laser 
linewidth. This quantity will be the focus for the remainder of this chap- 
ter. Linewidth measurements are carried out with the aid of a spectronv .er. 
Spectrometers come in a wide variety of types and operating principles. 

One type, the Fabry-Perot spectrometer, fits both the given wavelength 
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resolution and available equipment requirements. Section 3.2 will describe 
and present the theory behind the operation of a Fabry-Perot spectrometer. 
Section 3.3 will discuss the c .teria taken into consideration when 
designing an interferometer, while Section 3.4 reports the results of the 
linewidth investigation of the IJrbana lidar transmitter. 

3.2 Kil'ty-Perot Theory and Principles 

Fabry-Perot instruments operate on the principle of multiple beam 
interference. Multiple beam interference is the process in which an incident 
beam undergoes many divisions in amplitude, each division experiencing a 
certain amount of phase shift. The divided beams are then recombined 
allowing them to interfere with each other and form a characteristic fringe 
pattern. Figure 3.1 diagrams the interferometer used in this study. The 
two interferometer mirrors are the key members of the apparatus, since this 
is where the beam division is performed. Plane mirrors were used in this 
case, but there are other mirror configurations, con focal for instance, 
which also work. To simplify the discussion, only the plane Fabry-Perot is 
considered. 

Figure 3.2 shows a detailed view of the mirror surfaces and the ray 
geometry for one incident monochromatic beam. As the figure shows, a set 
of rays emerges from the mirrors, each successive ray having undergone a 
phase delay that will be called 6. The group of rays which forms the re- 
flected light from the interferometer won't be analyzed since only the 
transmitted rays make up the information that is being used to deduce the 
laser linewidth. Few oend Wolf [1959] provide a good analysis of these 
reflected rays. The delay, <5, is just the phase difference between ray 0 
and ray 1 as ray 0 propagates from points a to b while ray 1 propagates 
from points a to c to d. The delay for ray 1, 6^, is given by 
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where 4> is the phase shift encountered at each reflecting surface. The 
phase shift for ray 0, 6^ can be expressed as: 
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Combining these, 
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which reduces to: 


4 * cose 2$ 


(3.5) 


An expression for the transmitted electric field can now be developed. 
Assume the mirrors are identical, each having a reflection coefficient of r 
and a transmission coefficient of t. Using common phnsor notation, the 
electric field for ray 0 is given by: 


0 ■ V 


(3.6) 


where A () is the incident field strength amplitude and point b is taken as 
the zero phase reference. ' l's electric field can be written as: 
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(Generalizing the expression for ray n , wo have 

K 0 

The sum of all of these fields gives the total field expression: 
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where it is assumed that an infinite number of beam divisions take place, 
liquation (3.9) is a geometric series and sums to: 


‘V 

1 - r 2 . -«'* 


(3.10) 


Since detection devices for light measure the light’s intensity, the value 
for the transmitted intensity should be found: 
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If the incident intensity is denoted by and the reflection and 

t ransmission coefficient are replaced by the reflectivity, r ? * , and the 

transmissivity, (3.11) becomes 


‘0 


(3.12) 


1 + F" - 2K cos 6 

Iquation (3.12) is one of the forms of the Airy function. The transmitted 
intensity pattern will be dependent on the plate reflectivity H and the 


, * 
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phase delay <5. 6 itself is dependent on the beam's angle of incidence 0, 
the wavelength X , mirror separation d and phase shift at the mirror's sur- 
face #. Figure 3.3 shows a plot of the Airy function. 

A closer look at -5 is now in order, to see where constructive inter- 
ference will occur. Looking at equation (3.12), 7^ is clearly maximized 
when cos<5 is also a maximum, or when 6 is a multiple of 2r . That condition 
can be stated as follows: 

6 * -p- cosO + 2<p ■ 2qir (3.13) 

q being an integer. Next', incorporate the phase delay from the mirror 
reflection into q, making q generally a non-integer, that still increments 
by integers. This is done to simplify the ensuing manipulations and allows 
(3.13) to be rewritten as: 

(S r cosO ■ 2qn (3.14) 


q is usually a very large number for common interferometer configurations. 
Assuming a mirror separation of 1 cm, a wavelength of 500 nm and small 
angles where cosO can be approximated as 1, has a value of 4 x 10 . 
Therefore, very small changes in either cosQ, d or X will advance one to 
the (q + l)th order of interference. 

Figure 3.1 shows that the information recorded by the film in the 
interferometer is just the intensity variation as a function of 6. Figure 
3.4 shows a plot of the Airy function where only 0 is varying. This plot 
was made with i * 5 cm, F * 0.85 and X ■ 589 nm. The condition on 0 for 
intensity maxima is derived from equation (3.14) and is given as 

2d cosO » q\ (3.15) 

Note that q must decrease as 0 increases, which means that intensity maxima 
which occur at larger angles are lower orders of interference. To see how 
the angle of incidence for constructive interference depends on the order 
of interference, let q * - n where n is an integer ranging from 0 to 


<?0* 


infinity. Also, assiwie that -v-j- * 1 so that the "zeroth" order of inter- 
ference occurs along the optical axis, or at 0 * 0. Hquation (3.14) can 
then be written as 
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Figure 5.4 Plot of the transmission through a F ah r> -Perot \ersu^ angle. 

Note that the transmission peaks tend to group together for 
larger angles. 
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Since 0 will be small, cosfi can be replaced by the first two terms 

n n 

in its scries expansion, 


1 



nX 
* T7 


(3.17) 


Solving for 0 ^ then gives 


0 » (3. 18) 

n 1 

*) varies as the square root of the number of orders away from < 7 ^, This is 
confirmed by Figure 3.1 which shows intensity maxima grouping closer together 
us higher angles are reached. Note also that as the mirror spacing •! in- 
creases, 0 decreases indicating that the intensity variation pattern will 
he pushed together for larger mirror separations. 

The interference pattern that forms on the recording film will next be 
investigated. Figure 3,3 diagrams the hack half of the interferometer. 

The divided beams leave the mirrors at an angle that will allow for con- 
structive interference when recombined. The focusing lens provides the 
moons of recombination by focusing the rays to a point in the lens' focal 
plane. There will also he light focused at other points where a ray is at 
an angle with respect to the optical axis. Therefore, the total two 
dimensional pattern recorded by the film will be the point rotated a full 
revolution around the optical axis, or in other words a ring pattern. This 
is shown in Figure 3.6 which is an actual photograph of the ring pattern 
produced by the dve laser light passing through the interferometer shown in 
Figure 3.1. The radius of the ring for a particular 6 is easily found by 
looking at Figure 3.5 and is given by 


t an*' « -7 (3.19) 

or, for small " ; 

0 w *•; (3.20) 

Throughout the disci .ion, only a monochromatic source of radiation 
has been considered. However, the purpose of this study is to measure the 
linewidth of the incident radiation which implies a continuum of monochro- 
matic lines over a narrow wavelength region. Therefore, now consider two 
different wavelength components, one at the other a small measure of 





Figure 3.6 Photograph of the interference ring pattern of a laser 
pulse from the Urbana dye laser. The laser linewidth 
is proportional to the product of the ring's radius 
and thickness. 
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wavelength awnv at X - AX. If the same order of interference is considered 
for both, say q, then the angle of incidence for the X - AX ray must be dis- 
placed from the X ray to satisfy the condition for constructive interference, 
equation (3.141. Fquation (3.14) can be written for each wavelength giving: 
2 ./ COS A - q\ (3.21a) 

2d cos ( 0 ♦ AO) « < 7 (X - AX) (3.21b) 

where AO is the difference in angle needed to allow equation (3.14) to be 
satisfied at the different wavelength. Fxpanding the ens in equation 
(3.21b), it then becomes: 

2d cos 0 cosAO - Sine sinAO * q\ - qh\ (3.22) 

Several approximations are in order here. First both 9 and especially AG 
are small enough that their cosine may be approximated by 1 and their sine 
by the angle itself. This allows (3.21) to be rewritten as: 

■ q X (3.23a) 

2d - 0 M * q\ - qh\ (3.23b) 

Solving (3.23a) for q and substituting into (3.22b) then yields: 

= 6 AP (3.24) 


An expression for A0 can be found by differentiating (3.20). Substi- 
tuting this for AG and (3.20) itself for 0 gives an expression for AX in 
terms of measurable parameters r, Ar, f and X of: 


AX 


r Ar 


(3. 25) 


Ar represents the change in radius of the ring in going from X to X - AX. 

An expression has now been developed which relates the quantity under 
study, the linewidth AX, of the laser output to the information that's 
recorded by the interferometer, namely the ring pattern. Next, information 
on the design of an interferometer that will give meaningful results is 
needed. 

3.3 Fal'rit-rsj'ot Inter fevone ter Desian Criteria 

Re- exami nation of Figure 3.4 reveals that even for a monochromatic 
source, there is still a AG associated with the intensity peak. This AO 
would imply a corresponding A r which in turn would indicate the source would 
have a finite linewidth of AX as calculated by equation (3.24). This appar- 
ent AX is due to the non-unity reflectance of the mirror set and indicates 
a minimum resolvable linewidth of the interferometer system. Figure 3.3 
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shows the way the peaks narrow as the reflectance goes up, however, even 
perfectly reflecting mirrors would not allow for an infinitesimally small 
resolvable bandwidth. This is due to the amount of light transmitted by 
the interferometer going down as the reflectance goes up. Clearly in the 
limiting case of unity reflectance, there will be no light transmitted and 
hence no pattern to record and measure. As it turns out, there are factors 
other than the mirror reflectance which cause broadening in the intensity 
peaks. Therefore, one cannot get arbitrary resolution by choosing high 
reflectances and detectors sensitive enough to reconi the associated reduced 
intensity. The factor which plays an important role in determining the 
minimum resolvable linewidth is the mirror flatness. This topic will be 
discussed later, but first consider only the effects of the mirror 
reflectivity. 

Figure 3.7 shows a plot of the Airy function, equation (3.12), as a 
function of the phase difference < 5 . If only the wavelength is allowed to 
vary in the expression for 6 , which is given in equation (3.14), then the 
plot becomes a function of wavelength. As the plot shows, the Airy' function 
is quasi -periodic in wavelength, with a period that is called the free 
spectral range (FSR) of the interferometer. This is an important design 
parameter for interferometers and an expression will now be derived for it. 

First write the condition for peak intensity', equation (3.15) for the 
?th order and the (.7 - l)th order of interference: 

2 • coso = ,;\ (3.2<>a) 

2,; COSO = ( : - 1) (A ♦ AX pSR ) (3.2<>b) 

where AXp^ is the change in wavelength necessary to compensate for the 
change in order. AXj,,,., is the free spectral range of the interferometer. 

To simplify the derivation assume the ring pattern is observed at a small 
angle, allowing cosO to be replaced by one. Subtracting equation (3.26a) 
from (3.26b) gives, after some rearrangement: 

X . (, - 1) AX fsr (3.27) 

The small angle approximation allows equation (3.26a) to be solved for 
.7 giving 

*7 = 2. I/A (3.28) 

7 is going to be a very large number, as discussed earlier, allowing (<7 - 1 ) 
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to he replaced by i in equation (3.2 7 ). Incorporating both of these obser- 
vations into equation (3. 2"’) allows for the free spectral range to be solved 
in terms of measurable quantities: 


AX 


1'SR 


X“/2.. ; 


(5.29) 


Note that as the mirror separation is increased, the FSR becomes smaller. 

With this background, an expression for the minimum resolvable bandwidth 
of a system can be determined. Let A6 be the phase difference necessary to 
move from the maximum to the half maximum of the intensity plot. Writing 
the normalized Airy function at the half maximum point yields: 


1 (1 - /') 2 
T = “ 

1 + - 2r cos(5 0 + A A) 


( 3.501 


where A = 2j". Rearranging and cancelling common terms transforms (3.30) to 


-> 

cos ( A + '.61 = 4 - 1 - }'*“ 


(3.311 


l.xpanding cos(»S () + A A) and applying to the proper terms then gives: 


2’> cosAA = l : ' - 1 - 

Next, replace cos A A by the first two terms in its series expansion 



(3.32) 

(3.33) 


which reduces to: 


AS = (1 - ’■’) (3.34) 

v r 

Now, if AA is caused only by a change in wavelength, it can be expressed by 
differentiating conation (3.14) with respect to 1 giving: 


! A A ! = 


4 7i ’ cose 
X" 



(3.35) 


where AX, represents the change in wavelength corresponding to the change in 

phase A A . This will onlv be half of the minimum resolvable bandwidth AX . . 

nun 

Using this along with equations (3.34) and (3.35) allows the solution for 
the minimum resolvable bandwidth to be obtained: 


AX . 

min 


(i - ■■•■) x : 
v 2ni/ cos0 


(3.36) 


Making the small angle approximation and rearranging yields 
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(3.37) 


The numerator of this expression is just the FSR of the interferometer. 
The denominator is called the finnese, !', where 



For a Riven FSR, a larger finesse will Rive the interferometer a smaller 
minimum resolvable bandwidth. Therefore, the finesse can be regarded as a 
figure of merit for the F'abry-Perot device. Note that as the reflectance 
goes to unity, the finesse goes to infinity. 

This expression for the finesse was developed for perfectly flat 
mirrors with an arbitrary reflectivity. As mentioned before, a deviation 
from perfect flatness will cause the intensity peaks to broaden. The effect 
of this deviation from flatness can be expressed in terms of another finesse 
figure, the flatness finesse. Equation (3.38) then refers to the reflective 
finesse. For mirrors which deviate from parallelism by an amount \/M, the 
flatness finesse will be given by 

F f « M/2 (3.39) 

For mirrors flat to X/50, this would imply a flatness finesse of 



The total instrument finesse can now be calculated for real-life 
interferometers as 

(3.40) 

where F„, is the total instrument finesse and is the reflective finesse, 
equation (3.38). Therefore, the part of the finesse which is lowest will 
tend to dominate the total finesse. 

The interferometer design will be considered now. The two major design 
criteria are: 

1) To choose a free spectral range which is larger than the expected 
linewidth of the light to be measured, 

2) (liven the free spectral range from above, to choose a finesse large 
enough to give a minimum resolvable bandwidth smaller than the 
expected linewidth. 
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If condition one is not met, the rings will be so wide that they will over- 
lap orders and destroy the ring pattern. This is apparent after some 
thought about the meaning of the FSR. Condition two ensures that we are 
actually measuring the linewidth of the laser and not iust the instrument 
bandwidth . 

Kith an FSR chosen, using equation (3.29), the proper mirror spacing can 
bo chosen. Next, with the minimum required finesse Known, one can choose a 
set of mirrors that will have a large enough reflectance and sufficiently 
good flatness figure. In order for the mirrors to be aligned well enough 
to vield a high finesse, a mount that is very mechanically stable and is 
resistant to thermal fluctuations must be used. Commonly, some type of 
invar is used in this mount's construction to provide the thermal stability 
required. If the mount vibrates or expands over even small fractions of 
the light's wavelength, alignment suffers and the finesse is reduced. 

A review of Figure 3.1 also shows values of two lens focal lengths to 
be chosen. The focusing lens focal length should be chosen large enough to 
make the ring pattern large enough to easily observe, liquation (3.25) shows 
that for a given A 1 and ami A?* must increase to offset increases in 

thus providing a larger and more observable ring pattern. 'I he field lens 
focal length must he small enough to allow for a wide variety of incident 
angles so that plenty of orders of interference will be excited. If A is 
the laser beam diameter and the field lens focal length, it is easily 
shown from equation (5.18) that the number of rings excited by the field 
lens is given by 


Under ivmal conditions, •. is quite large and needn't be worried about, 
o.l .".’a ' n : , r. ' Von.’ V 

Previous evidence indicated that the linewidth for the IJrbnna dye 
laser would be somewhere around 10 pm |o •'.»»: fop :v •' , 19”8]. There- 

fore, using the design criteria discussed in the previous section, it was 
felt that an FSR of 15 pm would he appropriate for this study. This value 
for the FSR corresponds to a mirror separation of l.lo cm. The next step 
was to find a sot of mirrors which would provide a high enough finesse, 
when aligned properly, to provide the instrument with a low enough minimum 
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resolvable linewidth to perform the measurement. A set of 85 percent re- 
flective mirrors which were flat to X/20 were readily available. This set 
of mirrors has a reflective finesse of 19.3 and a flatness finesse of 10, 
yielding a total finesse, as calculated by equation (3.40), of 8.9. Using 
equation (3.3?) for this finesse figure and the previously mentioned FSR of 
15 pm, a minimum resolvable linewidth of 1.7 pm was calculated. This would 
work quite well for the estimated linewidth. These quantities and equipment 
specifications are listed in Table 3.1. 

It was important to have a good mechanically and thermally stable mount 
which would hold the mirrors and allow them to be adjusted to an accurate 
degree of parallelism. Figure 3.8 exhibits a photograph of the mount used 
in this study. The spacing rods are constructed of steel, with the mounting 
plates manufactured at close tolerances to reduce slipage on the rods. 

Again, Figure 3.1 diagrams the equipment setup. A description of the 
alignment procedure for this equipment now follows. Using the beam from 
the continuous wave HeNo alignment laser, the basic system equipment is put 
into place. This includes the iris, the folding mirror and the mount, 
without mirrors in place. The folding mirror is now adjusted so that the 
beam passes through the optical axis of the mirror mount. Next, the front 
mirror is put into place and it is adjusted so that the HeNe beam retro- 
reflects back through the iris. The hack mirror is then put into place. 

Some sort of dot pattern will most likely now appear on a viewing screen 
in back of the interferometer. This is the HeNe beam as it unstably walks 
off the mirrors. The back mirror is now adjusted to bring the dots closer 
together until they all fall on top of each other. This should be done 
with ns much precision as possible. The field lens is next put into place, 
two to three focal lengths away from the front mirror surface and positioned 
so that the beam stays in the optical axis. This can be done by placing an 
iris stopped all the way down in the HeNe beam at the point where the field 
lens focus is expected to fall. Then the field lens is inserted and ad- 
justed so that the focused spot falls in the center of the iris. This 
second iris can now be removed and the focusing lens put into place. Again, 
this lens should be placed so that its center falls onto the optical axis, 
although the adjustment is not so critical as for the field lens. A viewing 
screen is then placed one focal length away from the focusing lens. This 
distance should be measured with as much accuracy as possible to minimize 
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Table 3.1 Parameters and specifications of 
the Fabrv-Perot interferometer 
used to measure the linewidth of 
the Urban a dve laser. 


licjuipment/ Parameter 

Center Wave 1 en gt h (\) 

Reflectance (r’l 

P 1 at no s s 

Spacing ( :} 

Free Spectral Range 

Reflective Finesse (F. ,") 

Flatness Finesse (F v ) 

Total Finesse (F .) 

Minimum Bandwidth (A* . ) 

min 

Field Fens Focal Length (.•’„) 
Focusing Lens Focal Length (;") 


Value 

589 T.m 
.85 
y:o 

1.16 cm 
IS pm 
19.3 
10 
8.9 

1 . " pm 
S cm 
1(H) cm 
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t’uuj^iHr error. A ’■j'.ty pattern should, now In* observable . If the pattern 
appears unsyrimet ri c , ?.«•«.»» 1 adjustments can l»e made on the hack mirror in 
the interferometer set to a l low for the in- si. att ainable ring pattern. The 
covered recording film can now he put into place. The position of the film 
shou] 1 he in the focusing lens focal plane, i he dye laser Is then turned 
on .it this point to ensure the ring patten; has good symmetry and is posi- 
tioned proper 1\ with respect to the recording film. The apparatus, from 
the iu» er'Vr. 'meter mount on hack, should now he covered hy a light tight 
cloth to prevent stray light from reaching the film. The apparatus, after 
alignment, was photographed and is shown iri figure 5.9. 

To take a picture of the ring pattern, the beam to the interferometer 
is Mocked off until the film has been uncovered. One laser shot is allowed 
to oar.s through and then the film is re-covered and developed. This study 
simpiy used the hack end of an oscilloscope camera with type 40 Polaroid 
film, allowing quick, viewing of the photograph. If the photograph has not 
received tho correct exposure, the beam intensity must he adjusted. Tor 
over exposures , the use of a less reflective mirror in sampling the beam or 
an interference filter placed in the beam work well in reducing the incident 
energy on the film. The interference filter is especially suitable since 
its transmission characteristics are a function of the beams angle of inci- 
dence to the surface, Therefore, by varying the angle of the filter with 
respect to the beam, a continuum of intensities can be produced. 

figure 5." shows ;i photograph properly exposed, A trawling microscope 
was used to record the position of several of the rings for various cross 
sections through the ring pattern. Values for v times An were determined 
for the rings sampled, through each cross section. These were averaged 
together ami inserted into equation (3.2 5) along with f ■= 100 cm and * * 
nm. This vi elded a calculated linewidth of b.O * 2.5 pm. 
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SMOOTHING IN Till RAYL1: ICJI SCATTERING REGION 


4,1 •tioK 

In Chapter 2 , the following expression was developed Cor the average 
sodium concentration in a range bin '• 

Ha, . + Mi/ 2)" r ; , V. 

T ‘ *■ (A 


Note that the return from a bin in the Rayleigh scattering region, V, , 
appears in the expression. There are two criteria to Ik* considered when 
selecting an altitude at which to normalize the photocounts returned from 
the sodium layer, first, the altitude must be high enough that only a 
negligible amount of light is attenuated from the beam as it propagates from 
the normalizing altitude to the bottom of the sodium layer. Secondly, the 
returns from this region must be due to Rayleigh scattering. This is neces- 
sary since Rayleigh scattering is ordinarily the last loss mechanism operating 
on the be;im until it reaches the sodium layer. Also, Rayleigh scattering is 
well understood and expressions which relate the returns via that mechanism 
are easily obtained. 

This chapter will deal with issues related to choosing, if possible, a 
range and an associated bin value to normalize the returns and calibrate the 
photocounts. Section 1.2 discusses the predicted nature of the Rayleigh 
returns. Section 4.o mentions problems due to photomultiplier tube (PNfT) 
overloading from this region. Section 4.4 presents a method to smooth the 
selected range bin to reduce errors contributed by noise. 

4.2 ;•> f U •*>:«■» *y>rr iiw /•.:;< ’rf-ih 

\n expression for the returns from an altitude where Rayleigh scattering 
is the dominant scattering mechanism was developed in Chapter 2 and is given 
by: 

*> ^ 

•*VVV“ r » h;. *'•" n t: .(-n) 

V., u — ^ (4.2) 

The terms V, , .4 ( , n and ,\r. will be constant with height. Two terms, ? K and 
j (tt) vary slowly with height and can be considered constant over the inter- 
val of altitudes under scrutiny here. The average atmospheric density term 
~ follows closely to the following exponential relationship: 
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V* 


'R °o* " C4.3) 

where is the density at an altitude of zero and H is the atmospheric 
scale height. Lumping the terms which don't vary with height into a constant 
K, equation (3..?) can be rewritten as: 


Nr 


K 


R 


-V* 


(4.0) 


Over a narrow altitude range the inverse altitude squared relationship 
will be dominated by the exponential. Therefore, returns from the Rayleigh 
region should demonstrate a strong exponential decay with height. Figure 
4.1 demonstrates returns from the Rayleigh region showing this predicted 
exponential response. 

4.3 Photomultiplier Tube Overloading 

With the addition of the 1.22-m telescope to the Urbana Lidar system, an 
increase of a factor of 5 over the 0.38-m telescope has been observed in the 
photocounts from the Rayleigh region [Rowlett md Gardner, 1979]. This in- 
crease in signal has been significant enough to cause the system PMTT to over- 
load. The overloading is characterized by returns much reduced or even nulled 
out over a period of time followed by a recovery interval where the counts 
rise to their expected value. The null or reduced response period depends 
on the strength of the overloading signal. Figure 4.2 depicts a plot of the 
unfiltered photocounts vs. altitude of a profile taken with the 1.22-m tele- 
scope. The returns show a null portion up to the 25 km range followed by 
the recovery interval up to approximately 40 km after which the returns 
follow their usual response. Clearly, there is no exponential decay observ- 
able in the lower regions which means that this profile cannot be calibrated. 

Fortunately, the overloading is not always so severe. Figure 4.3 shows 
another profile of unfiltered photocounts plotted vs. altitude. The region 
below 30 km shows erratic anJ reduced returns until reaching the 28 km level 
where an exponential decay is clearly evident. By selecting a bin in the 28 
km to 45 km interval where Rayleigh scattering is dominating the returns, 
the profile can be calibrated. The following section presents a method for 
choosing this altitude. 

4.4 Rayleigh Normalization Bin Selection and Smoothing Technique 

Figure 4.4 is a plot of the photocounts vs. altitude for a profile 
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which has an interval strongly dominated by Rayleigh scattering. Clearly 
some amount of noise is present in this signal as the deviation from a 
purely exponential decay reveals. To minimize the effects of this noise, 
some sort of smoothing needs to be done on the Rayleigh scattered bin 
values. The following presents a least squares fit method of performing 
the desired smoothing. 

Equation (4.4) expresses the number of photocounts received at a height 

n 

r . v as a function of that height. By multiplying each side of (4.4) by 
and taking the natural logarithm of each side, (4.4) now becomes: 


In 



In 


A " H 


(4.5) 


liquation (4.5) exhibits a linear relationship between In n f z r and the alti- 

tude a p . Therefore, a least squares fit can be done for these two quantities 
providing the desired smoothing. 

As discussed previously, not all of the returns are due to Rayleigh 
scatter. Range bin values which are dominated by noise or which show signs 
of overloading effects in the PMT cannot be used in this fit. A computerized 
method to examine the returns and choose and smooth an appropriate bin is 
used in this calibration technique. In this method a background noise level 
is first computed by averaging the returns in several bins which are due to 
noise. Two thresholds are now computed. An upper threshold at four times 
the background level and a lower threshold at two times the background level. 
The returns within these threshold levels are examined and the altitudes 
where the returns are at a level of three times the background counts are 
found. All pertinent levels and thresholds have been indicated for the 
profile in Figure 4.4. The altitudes at which the count is three times the 
background count are averaged to determine the height at which the normali- 
zation is performed. Next, the five bins on either side of the chosen bin 
are curve fitted and a smoothed value for the normalization bin is computed. 
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5. liXPHRIMFNTAL RESULTS AND HRROR ANALYSIS 
5 . 1 Tx/v ft'wn ta.1 He an 1 1 a 

In Chapter 2, an expression relating the photocounts to the sodium con- 
centration was developed, the final form of which is given by equation (2.45). 
Chapter 5 provided the necessary information to calculate a value for the 
effective cross section. The value of the Rnyleigh range bin, its altitude, 
and the number density of the atmosphere were the topics for Chapter 4. 

With the addition of the raw data which provides values of altitude and 
range bin photocounts, equation (2,45) can be evaluated. These results are 
the topic of this section. 

As mentioned in Chapter 1, a seasonal variation in sodium column abun- 
dance is defined as the number of sodium atoms contained in a column extending 
from the ground to the top of the atmosphere and having a unit area cross 
section. The column abundance reported in the past shows an increase in the 
winter values over the summer values. This feature is evidenced by comparing 
figures 5,1 and 5.2. Figure 5.1 is a calibrated profile from the June 24-25, 

1979 data set. Note that the layer extends from 80 to 100 km with a peak 

0 - 5 

concentration of 1.5 * 10* m * at an altitude of 85 km. figure 5.2, the 

winter profile was taken January 18-10, 1080 and extends over the same 80 to 

100 km range. The peak of the layer now occurs around 88 km and has a value 
0-5 

of 7.5 " 10* m *. Most of the increase in sodium atoms is found on the top- 
side of the layer. libit in ! ;■ [19T] have also observed this 

occurrence . 

Figure 5.5 illustrates the average nightly column abundance against the 
month of the year. The data was taken during a period extending from 
November, 19'’" to January, 1080, examination of this figure shows some 
interesting features, first, the abundance in the winter months of December, 
January and February are a factor of 5 to 7 times larger than the abundances 
observed during the summer months of June and Ju’y. Note also that the 
cyclic nature of the variation is fairly smooth. In other words, approxi- 
mately equal amounts of time are spent in the enhanced and depleted portions 
of the cycle. This smooth variation suggested that the data lie fitted, by 
means of a least squares error method, to a sinusoid. The functional form 

of the sinusoid is given as 5.75 x 10**' sin(27T(x - 9.85)/12) + 4.55 * 10**^ 

*> 

m~ ~ where x is the month. The value of x ranges from 0 to 12 with 0 


SODIUM CONCENTRATION (IO»rrf s ) 

A calibrated profile of summertime data taken on the 
evening of dune 24 and 25, l<r<>. This particular profi 1 
was t.ie fifth one taken that evening and was low-pass 
filtered with a cutoff spatial frequency of <5.~5 kmr 1 




ALTITU 




OBSERVATIONS WITH 5 OR 
LESS PROFILES 



i lot showing the seasonal average o' the ^odiun column abundance, 
data were fitted via a least-square error method to a sinusoid havi 
a value of 5.” \ M 1 -' sin|:-|.V - 9.85; 12| * 4.55 \ If) 1 --' n'-. ,he 
tern .V is a variable representing the month. \ value of 0 represen 
the beginning ot .January while jne ot 15, the end of I'ecen'* , cr. 






representing the beginning of January and II the end of Oeeomher, As can be 
seen the data fits the curve well although a larger data base would be needed 
to draw more definitive conclusions, The phase delay of 9.85 corresponds to 
October J(». 

While these results compare well to the past observations reported by 
other groups on the whole, there are some differences. First, the factor of 
5 to ” increase agrees much more closely with the th'l'jcK a>ui dan i (1971 j 
value of 5 than it does with the AWz\ ( 1 7 7 ] value of 4 and the 

# [I9"9j value of J, d . [|9~9| felt that this 

factor was related to latitude, larger factors being found at higher lati- 
tudes. The latitude of the Wrbana system of 40°N and the Menie dUm^nt 
(197” | latitude of 44°N are in close agreement. However, the factor of 
increase agrees more closely with the fnglish results obtained at a latitude 
of f.2°\. Therefore, the latitude dependence of the factor of increase does 
not seem to he valid for these results. 

.W r* d. [19”9J reported a smooth cyclic variation, however, 

AV|,\V v: ; id ;rx (1977) June observed a much shorter period of enhancement 
lasting only 2 to S months, the rest of theyear being somewhat equivalent 

jk : .%'nw [1971| published a variation similar to this one but having 
a marginally shorter enhancement. 

finally, the consistency of results as shown in Figure 5.3 yields con- 
fidence in the relationship between the calibration of the different evenings 
of observations, So points, even those with less than 5 profiles averaged, 
deviate significant 1\ from the rest of the group. 

5.J f ro, i' 

io perform the error analysis, an RMS mean error will be used, This 
can be defined as* 




( 5 . 1 ) 


where a delta quantity represents the uncertainty of that quantity. Assuming 
: , r,., \' v , , and e, (*) have negligible amounts of error, p^ becomes 

a function of , _ , A. and a <v ,. Ivaluating equation (S.l) using these varia- 
bles results in: 
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relation of the profiles to each other should be much more reliable, However, 
for this tuning error, the relationship between profiles becomes unsure as the 
tuning randomly varies through the evening. 
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Table 5.1 Summary of statistics from the October 29, 1978 dataset. 

The column abundance is directly related to the signal 
photons and inversely related to the Rayleigh photons. 
Low numbers of signal photons m profiles 4 and 6 are 
probably due to drift in the laser tuning since the 
Rayleigh counts stay relatively constant. 


Profile 

# 

Signal 

Photons 

Rayleigh 

Photons 

Column 

Abundance 

\ Deviation 
from 

avg. abundnn 

1 

JO 8b 

30.4 

6.3 

X 

io 13 

26 

2 

171b 

35.6 

6.3 

X 

10 13 

2b 

3 

1263 

31.6 

4.6 

X 

io 13 

- 8 

4 

730 

28.2 

3.3 

X 

10 13 

-34 

5 

144b 

31.2 

5.4 

V 

10 13 

8 

6 

872 

33.5 

3.7 

X 

io 13 

-26 

n 

t 

105b 

24.8 

H . 3 

X 

10 13 

-14 

8 

1240 

33. 7 

5.1 

X 

io 13 


9 

1225 

30.0 

6.1 

X 

10 13 

22 

10 

1551 

28.0 

6.9 

X 

10 13 

38 


Avg 
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6. SUMMARY AND RECOMMENDATIONS FOR FURTHER WORK 

6 . 1 Fi mnary 

From the lidar equation, an expression was developed to calculate the 
sodium concentration from observed photocounts collected with the llrbana 
lidar system. The main unknown in the equation was the laser linewidth. 

A Fabry- Perot interferometer was designed and constructed to perform this 
measurement. The laser linewidth was measured to be 6 t 2.5 pm. Some 
problems in calibrating the photocount vs. altitude data were encountered 
from photomultiplier tube overloading. A method for selecting a normaliza- 
tion bin found in the Rayleigh scattering region was developed. With this 
information, past data were calibrated. The data showed a strong winter 
enhancement in the sodium column abundance of a factor of 5 to 7 over the 
summer values. The data were fitted with a cosine curve with good results. 
The calibration is accurate to within 50", approximately. 

6.2 Fe commend at forte for Further Work 

6.2.1 Further lineviith studies. A more rigorous study of the line- 
width and its variability needs to be undertaken to reduce the uncertainty 
in the calibration. Variability studies through an evenings* run and on a 
night -to-ni ght basis would be beneficial in ensuring an accurate calibration. 
An even better arrangement would be to monitor continuously the linewidth 
and record its values throughout the evening. Some sort of linear detector 
array would be best since it must have a quick response time. If photo- 
graphic methods are continued, some sort of alternative to Polaroid film 
should be used. The film should be able to accept a wide range of pattern 
intensities to cope with the changing output energies. 

6.2.2 Center frequency monitor. As was discussed in Chapter 5, 
drifting of the laser center frequency is another cause of error which 
should be eliminated. This parameter will need to be controlled and 
monitored during data collection. Either an interferometry or cell absorp- 
tion method shoulu be used. 

6.2.5 Gatina or shutterina of photomultiplier tube. In order to 
eliminate the PMT overloading, the PMT must either be gain switched elec- 
tronically or have its input mechanically blocked by a shutter during the 
low altitude returns interval. Electronic gain switching will require 
modification of the PMT housing as well as the external switching and timing 


75 


circuitry. A rotating shutter would require a motor equipped with a 
propeller-like shutter as well as position sensing and timing circuitry. 

This is a significant problem because some photocount vs. altitude profiles 
cannot be calibrated due to the severity of the overloading. 

0.2.4 Larger data base of calibrated \:ofiles. In order to achieve 
more confidence in the seasonal fluctuation of sodium, a larger data base is 
needed, especially from January to June. This information can be had at no 
loss to other ir vest i gat ive objectives. 
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APPENDIX I 

Operating Procedures for the Lidar Transmitter 
The following is a step-by-step instruction set for aligning and 

operating the Urbana dye laser. 

1. Prepare and install a new set of flash lamps. 

2. Prepare a concentrated solution of the Rhodamine 6G dye. 

7 >. Using the Simer pump on the stand holding the water deionization 

system, pump the cooling water from its reservoir into the nalgene 
holding tank. Using acetone, clean out the cooling water rvoir 
and paper towel dry. Be careful not to inhale the acetom . umes . 

Replace the cooling water reservoir and open up the stopcock at the 
bottom of the holding tank. Keeping the drain valve closed on the 
deionizer, open the inlet valve to obtain a good water flow through 
the deionizer. Slower flows haven't proven to increase the waters 
resist ivity . 

4. Remove the fire bricks which form the oven for the sodium resonance 

cell. Remove the cell and look at the base to see if there is an ade- 
quate supply of sodium deposited on the glass. If there is a sufficient 
amount, one should observe a silver-blue coating around the glass that 
obscures the sodium channels within. If the coating is not very thick 
and allows easy viewing of the sodium channels, more sodium needs to be 
released from the channels. To accomplish this, obtain a low voltage, 
high current ac power supply capable of an output of 9 to 10 amps. 

Notice that the cell has three wires exiting the bottom of the cell. 

Each of the outer wires corresponds to one of the sodium channels. The 
center wire is common to both. When beginning with a new cell, use 
only one of the channels to release sodium until it has been used up, 
then switch to the second. It is therefore advantageous to mark the 
lead of the channel which has been used most recently. Connect the 
power supply to the center lead and the proper outer lead and turn the 
power on to the supply. Begin at zero current and increase slowly to 
nine amps. Keep the current on until there is a good deposit of sodium 
on the glass. This should take 20 T: 10 seconds depending on the 
amount of use the channel has experienced. Replace the cell in the 
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oven and turn the heater on to a setting of 55 on the variac. Monitor 
the temperature closely, not allowing temperatures greatei than 150 C. 
Temperatures greater than this c^use excessive decay of the ceil. If 
time is short the variac can be turned up higher to get the oven to 150 
0 more quickly, but should then he turned back down again towards the 
55 setting. The cell should be positioned inside the oven in such a 
way that will allow the laser beam to enter through the oven window, 
pass through the center of the cell and out the oven's back window. 

Turn on the sodium arc lamp and adjust both it and its associated 
focusing lens to image the pinhole onto the entrance slit of the spec- 
trometer. Check to see if the sodium doublet is well defined and 
bright. The eyepiece of the spectrometer will need to be adjusted to 
form well defined lines. If the sodium doublet is not bright enough to 
see clearly, readjust the focusing lens and arc lamp. 

Drain any old dye which may be in the dye reservoir. This can be done 
by closing the valve on the side of the reservoir rack which leads back 
into the dye flow system. Afterwards, open the other valve which leads 
to the free drainage hose. Insert the drainage hose into a container 
capable of holding the amount of dye in the reservoir. On the front 
panel, flip the switch under the amber light to the on position. The 
dye should now he flowing into the drainage container. After all the 
dye has been removed, turn off the switch and rinse the reservoir with 
ethanol to remove any residue. If would be wise to wear gloves so that 
contact with the Rhodamine bG dye will be minimized. Turn vhe switch 
hack on again and drain out the ethanol which was used to rinse the 
reservoir. .Shut off the drain valve and replace the drain hose to its 
storage position. Reopen the valve which leads hack into the dye flow 
system. The dye filter should now be emptied. There are two ways to 
do this. First, if the dye filter is known to he fairly clean, place 
a container of at least 600 ml volume beneath the dye filter and open 
the wing valve at the bottom of the unit. Turn the pump to the on 
position which should cause the dye to empty out. When finished, turn 
off the pump and close the wing nut valve. If the filter is dirty ami 
needs replacement, again place the container beneath the dye filter 
and open the wing valve. Using the wrench kept on the side of the 
reservoir rack, loosen the bolt on the top of the filter. As the dye 
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filter container separates from its holder, dye will begin to flow into 
the vessel which is holding the discharged dye. Allow enough dye to 
drain out to allow the dye filter container to be comfortably threaded 
out of the dye flow system apparatus. The wound filter can now be 
accessed and replaced if necessary. To empty the dye in the dye rod 
and the tubing which connects it to the pump, locate the hose which 
enters the front panel of the dye flow system. Placing a beaker or 
other suitable container beneath this hose, remove the hose from its 
fixture and allow the fluid in it to gravity drain out. 

7. Refill the dye reservoir with new solvent. The mixture ratio is 1 part 
water to 1 part ethanol. Normally 95% ethanol is used, however purer 
percentages of ethanol can be used. The normal amount of solvent used 
for a data run is 16 liters. This would correspond to 8 liters of 
water and 8 liters of ethanol. When the solvent has been added, the 
pump can then be turned on. This allows for the air bubbles flapped in 
the system to work their way out of the flow system. To turn the pump 
on, flip the switch under the amber light on the pump control panel to 
the on position. The pressure should begin to build as will be evi- 
denced by the front panel gauge. When the red light goes out, the 
pressure has exceeded the threshold point at which the safety cutoff 
will shut down the flow system. Now return the switch under the amber 
light to the off position. The amber light should remain or. and the 
dye should begin flowing through the system. Next, add in the concen- 
trated dye mixture at a rate of 30 ml dye mix to 1 liter of solvent. A 
16 liter solvent solution corresponds to 480 ml of concentrated dye 
mix. Finally remove the 1, 3, 5, 7, cyclooctatetraene (COT) from the 
freezer and add it to the solvent at the rate of 0.2 cc per liter of 
solvent. Therefore, for a 16 liter solution, 3.2 cc of COT would be 
added. This is best done with a syringe and hypodermic needle. Gently 
stir the whole solution with a stirring rod. It will pr oably take 5 
to 10 minutes before the solution has mixed well with the dye flowing. 

8. Remove the f 1 ask that serves as a trap in the vacuum system. Fmpty any 
ablation product which might be in it. Periodically rinse it with 
acetone or ethanol and allow it to dry. Using a paper towel, remove 
the old vacuum grease around the flask's mouth and from its stopper. 
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Replace it with a fresh thin coat of vacuum grease on the contacting sur- 
faces. When the vacuum system is first turned on, this stopper may need 
to be held in place until the vacuum begins to suck it down on its own. 

9. If there is enough water in the cooling water reservoir to start the 
pump, now is the time to do so. Place the hoses in the water and after 
checking to make sure all hoses around the flashlamps have been recon- 
nected and the clamps on them removed, turn on the pump. Tf the pump 
does not prime, place a rubber bulb in the water and fill it with water. 
Keeping both it and the inlet hose under water, force water into the 
inlet hose until a solid column of water reaches and goes through the 
pump. The pump should now work at this time. 

10. Remove the cover for the bi-elliptical cavity and check to make sure the 
dye rod is oriented properly with respect to the cavity. It should be 
positioned along the common focus axis. This provides a rough alignment 
of the dye rou along the optical axis. 

11. Remove the optics from the rail and turn on the lleNe alignment laser. 
Adjust the laser so the beam passes through the center of the quart" 
windows on e: eh side of the rod. The dye should he flowing at this 
point since the beam propagates through the rod differently when there 
is no flow. An adequate way to center the beam in the windows now fol- 
lows. First, remove the pinhole from in front of the alignment laser. 
This is done to allow the lleNe laser to he moved without having to read- 
just the pinhole. Next, move the front end of the lleNe laser (the end 
nearest the cavity) to center the sopt in the window horizontally. If 
the vertical adjustment is off, check to see if the beam is at the same 
height above the rail for the length of the rail. Tf that adjustment is 
not correct, level the beam using computer cards underneath the alignment 
laser. When this is done, recenter the beam horizontally by moving the 
front end of the laser, as before. The beam should not be too far from 
the center of the nearest window. Undo the set screw which controls the 
height positioning of the rod and adjust the rod to center the beam in 
the window. Retighten the set screw and perform the same adjustment or. 
the hack window. The horizontal adjustment on the hack window will 

most likely be off at this point, but the window can still be verti- 
cally centered. Now, note to which side of the hack window the beam is 
positioned. Move the back end of the laser to the same side. The 
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amount moved will depend on the distance from center. Adjust the front 
end of the laser to center the beam on the front window. Check the 
position of the beam on the back window and make more adjustments as 
necessary. After completing the centering process, the bi-elliptical 
cavity may need to be readjusted so that the dye rod is still at the 
common focus and not in physical contact with the cavity. This is done 
by loosening the four bolts which hold the cavity to its black support 
frame, centering the cavity with the dye rod and retightening the bolts. 
Ensure that once the dye rod is centered, the hoses facing the dye rod 
will not be jostled or moved around. Even slight movements can destroy 
the alignment. This can be seen by moving the hoses and watching how 
the beam moves on a viewing screen at the back of the optical rail. 

12. Place the grating as far back as it will go on the optical rail. If 

it's not known whether or not the grooves are oriented vertically, they 
should be checked for proper alignment. To do this, first position the 
mount so that the beam is retroreflected back into the laser. Now, 
rotate the grating so that the angle of incidence 4 s changed. Soon, 
the first order of diffraction should appear somewhere around the dye 
cell. If the grooves are vertical, the first order too will retrore- 
flect back into the laser. If not, the grating should be rotated in 
its mount until this condition is met. Now turn the mount so that the 
first order on the other side cf normal incidence comes into view. 

This diffracted beam should also retroreflect if the grooves are prop- 
erly aligned. Make adjustments as required so that all three beams 
will retroreflect when the mount is rotated on the optical rail. Next, 
carefully examine the grating and choose an area, at least 5 mm in 
diameter, that is fairly clean and unmarked. Position the mount so 
that the laser beam is incident on this spot. Orient the mount’s angle 
to this beam so that the first order of diffraction is retroreflected 
back into the laser. The angular positioning micrometer should be posi- 
tioned at a setting of at least 300 mils, since this will need to be 
readjusted later. The mount should be clamped securely into place, and 
fine adjustments made with the micrometers to re-center the beam on the 
alignment laser. The pinhole should be in position in front of the 
laser to make the adjustment as accurate as possible. 
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13. Place the mirror that makes up the mirror grating combination in front 
of the grating. Select a good clean position on this and have the beam 
strike the mirror there. Securely tighten it to the rail and adjust 
the micrometers so that the beam is ret rore fleeted. This adjustment is 
critical and should be done with care. The grating may need to be 
blocked off so that its reflections will not hinder the mirror align- 
ment. 

14. Put the output mirror in place. Clamp it to the rail and adjust the 
micrometers so that it ret rore fleets. Usually, some faint spots, 
originating from the mirror grating combination, will be visible and 
can aid in this process. 

15. Between the output mirror and alignment laser, place a glass slide and 
position it so that the HeNe beam is reflected towards the beam split- 
ting table located on the sides of one of the 40 MHz transmitter cavi- 
ties. Position the beamsplitter mirror so that it guides the beam to 
the input slit of the spectrometer. The throughput from the beam- 
splitter should be guided via the total reflector through the center of 
the sodium resonance cell. 

16. Subtract 180 mils from the angular adjusting micrometer head on the 
grating. This converts the alignment from 632.8 nm to 589.0 nm. 

17. Block off the alignment laser. 

18. Turn on the water spigot on the mezzanine which cools the dye. Just 
a good steady flow of water should be sufficient. 

19. Replace both the bi-elliptical cavity cover and the wooden cover which 
goes over the whole central assembly. 

20. Turn on the vacuum pump, making sure the manometer's stopcock is tum.d 
so that both sides of the column are at the same pressure. As mentioned 
before, the stopper in the flask trap may need to be held down until the 
vacuum is well started. Allow this to run for about a minute so that a 
good vacuum is developed. Once a good vacuum is present, turn the stop- 
cock on the manometer to isolate the two sides of the column. This pro- 
vides the manometer with a vacuum reference to measure system pressure 
with . 

21. Turn on the nitrogen supply at the tank and adjust the diaphragm so 
that the line pressure gauge just barely shows some pressure. Open up 
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the line to the flashlamp system and allow the manometer to come to a 
reading of four torr. This valve should only be handled by the black 
plastic to avoid shocks when the high voltage is on. The shocks are 
not dangerous but do provide some discomfort. 

Turn on the pulsing equipment sitting on the table nearest the output 
end of the laser. There are two units which perform the pulsing. The 
high voltage pulscr is a metal chassis with a transformer mounted on top. 
The power switch and red indicator light are located on top. The second 
unit is aa external trigger unit in a small metal chassis box. Tins has 
the time between pulse adjustment on it, which should he set on the two 
second position. Turn the black slide switch on the high voltage pulser 
to the external trigger position and gaze into the bielliptical cavity 
to try to see if the flashlnmps are breaking down. If you cannot see the 
pulsing, listen for a snapping sound which will also indicate whether 
or not the pulscr is working. 

The high voltage is now ready to be applied. Quickly recheck to see if 
the dve is flowing, the cooling water is on and flowing, the dye cooling 
water is on, and that the vacuum system is on with pressure at four 
torr. Plug in the power supply and switch it on with the variac in the 
"zero voltage" position. Slowly increase the voltage until the lamps 
begin breaking down. This should occur around 13 kV. Continue the 
voltage increase until lasing is observed. The easiest place to look 
for lasing is on the grating. 

Adjust the up and down positioning micrometer on the grating to try to 
obtain a good symmetrical lasing spot. As you near the optimum setting 
the energy per pulse will increase. The voltage may need to be turned 
down to keep the grating safe from danger as the energy increases. When 
the spot resembles a low order spatial mode, gaze into the spectrometer 
and adjust the angular micrometer to tune the emission onto the line. 
The line is the right side of the sodium doublet looking into the 
spectrometer. Now, reduce power until the system is barely lasing. 

Make final adjustments on the up and down setting of the grating to 
optimize the spot configuration. If a ring or doughnut spot occurs, the 
laser is in good alignment and will probably work well. Don't use too 
much time trying to obtain this pattern, since it doesn't always occur. 
The laser will also work well if a good f.aussian spot is obtained. 
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25. Turn off the pulscr and reduce the high voltage power supply to zero 

and turn it off as well. Place the resistive discharge probe on top of 

the resistor bank and discharge the capacitor and power supply. You 

can watch the voltage reading on the power supply to see when the vol- 
tage has decayed to zero. Place the dead short on top of the resistor 
to get rid of any residual charge. Unblock the alignment laser and 
remove the wooden cover. Put the etalon in place and clamp it securely. 
As always look first to be sure a good clean area of the optics is 
being utilized. Adjust the etalon so that the laser beam is reflected 
to a spot just between the dye rod anil the flashlamp adapter on the 
bakelite plate. The beam should be at the same height as the center of 

the window. There should be a spot on the bakelite to mark this 

posit ion . 

2<i. Block off the lloN’c beam once more and replace the wooden cover. In- 
crease the system pressure to " or 8 torr. Switch on the pulscr ami 
look or listen for evidence of pulsing. Turn on the power supply and 
increase the voltage until good lasing is observed. The need for 
caution concerning damaging optical energy too high energy densities is 
reduced with the etalon in the system. 1'he optics have proven to be 
safe even at maximum voltage. 

2 ". book into the spectrometer and observe the laser emissions' spectral 
qualities. There should only be one thin narrow line somewhere close 
to the : line. Sometimes a faint line will occur one or two free 

spectral ranges of the etalon away from the main line. This is due to 
lasing off of the mirror in the mirror grating combination. Not much 
can be done with this problem. If other emissions are present, return 
to the grating alignment stage and begin again from there. Now look 
through the spectrometer and see if the laser line is close to the 
line. Under normal circumstances, it should look to be right on the 
line. Haze into the resonance cell and see if scattering is observed 
in the cell. If not, note the current position of the angulai adjust- 
ment micrometer and begin moving, the micrometer in stens of about 1 
mil in one direction, looking for scattering after each new adjustment. 
After moving 10 mils away from the center position, return to the 
original position, just noted, and go in the other direction in 1 mil 
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steps. Again, if no scattering is observed after moving «i0 mils in 
this direction, return to the 10 mil position on the other side of the 
starting point and search through 10 more mils in this direction. If 
no scattering has yet been found, look back into the spectrometer and 
adjust the micrometer to recenter the laser line on the line. This 
is the procedure to follow if the laser line was originally displaced 
from the line. Note this new position and restart the search proc- 

Am 

ess. If again no scattering is found, return to step 25 and repeat 
those procedures, perhaps trying a new spot on the etalon. Once you 
have found scattering, locate the positions where scattering is just 
barely observed. Make the final setting to be the middle of these two 
positions. Usually, there is a spread of 1 to 2 mils in these posi- 
t ions , 

28. Remove the glass slide from the rail and place the 45 degree mirror 

into place, clamping it down securely. Adjust the micrometers on this 
mirror to center the hoam on the marks on the hatch. Open the hatch at 
this point being careful of the high voltage. You may want to discharge 
the system if so desired. The laser is now ready to be aligned with 
the telescope. 

2‘>. As the tele. scope/ laser alignment is being performed, place a silvered 
mirror behind the 45 degree mirror and direct this residual beam 
through the sodium cell and into the spectrometer as before. Recheck 
the tuning every It) to 20 minutes and re center the micrometer between 
the "just scattering" positions. If pulsing becomes irregular, adjust 
the pressure upwards if spontaneous shots are occurring and downwards 
if expected shots don't occur. 

50. To shut down the system, turn off the pul so r, turn down the high vol- 
tage power supply to zero and switch it off. Unplug the power supply, 
then discharge the capacitor and power supply as mentioned before. 

Turn off the nitrogen and allow the column of fluid to come to equili- 
brium. Turn the stopcock to the horizontal position so that both sides 
of the column are at the same pressure. Turn off the vacuum system and 
then the nitrogen tank. Remove the nitrogen hose from the valve which 
is used to adjust the residual pressure in the vacuum system. Open it 
slightly and allow the flashlamps to return slowly to atmospheric 
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pressure. This should reduce the probability of leaks occurring in the 
flashlamps. Cover the optics with plastic bags, being careful the 
plastic does not touch the optics surface. Turn off the sodium cell 
heater and arc lamp supply. Turn off the cooling water to the dye. 
Close the hatch and tie it down. Turn off the dye pump by pressing the 
black button below the red lamp until the amber lamp goes off, then 
release it. Turn off the flashlamp cooling water pump. Finally, turn 
off the bench power strips. 
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